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ABSTRACT 

In the present article a new class of exact solutions of Einstein’s field equations for 

charged anisotropic distribution is obtained on the background of pseudo-spheroidal 

i+k4 

spacetime characterized by the metric potential g rr = -f-, where K and R are geo- 

l+ F 

metric parameters of the spacetime. The radial pressure p r and electric field intensity 

1. . ft - “jl n. Cl(K— 1)4 

E are taken in the form 8 np r = 44--— 2 —— and E 2 = —--The bounds of 

(■+**) 

geometric parameter K and the parameter a appearing in the expression of E 2 are ob¬ 
tained by imposing the requirements for a physically acceptable model. It is found that 
the model is in good agreement with the observational data of number of compact stars 
like 4U 1820-30, PSR J1903+327, 4U 1608-52, Vela X-l, PSR J1614-2230, Cen X-3 
given by Gangopadhyay et al [Gangopadhyay T., Ray S., Li X-D., Dey J. and Dey M., 
Mon. Not. R. Astron. Soc. 431 (2013) 3216]. When a — 0, the model reduces to the 
uncharged anisotropic distribution given by Ratanpal et al. [Ratanpal B. S., Thomas 
V. O. and Pandya D. M., arXiv: 1506.08512 [gr-qc](2015)] 

Subject headings: General relativity; Exact solutions; Anisotropy; Relativistic compact stars; 
Charged distribution 
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1. Introduction 


Mathematical model for generating superdense compact star models compatible with 
observational data has got wide attention among researchers. A number of papers have 
been appeared in lite rature in the recent past alo ng this direction considering matter distributio n 


incorporating charge (Maurya and Gupta 201 la,b 


c; 


Pant and Maurva 


2012 


has been suggested, as a result of theoretical investigations of Ruderman (1972) and 


M 


aurva et a/l20151). I t 


Canuto (119741) . 


that matter may not be isotropic in high density regime of 10 15 gm/cm 5 . Hence it is 


to construct charged distribution incorporating anisotropy in pressure. 


Bonner (1960 


pertinent 


19651) has 


shown that a spherical distribution of matter can retain its equilibrium by counter balancing the 
gravitational force of attraction by Coulombian force of repulsion due to the presence of charge. 


It was shown by 


Stettner ( 1973!) that a spherical distribution of uniform density accompanied 


by charge is more stable than distribution without 


on sp. 


te roidal spacetimes have been carried out 


1998), 


Sharma et al (2001), 


Gupta and Kumar (2 005 ), 


bi 


charge. The stud 1 

/of ch 

arge distributions 

Pate 

and Kopper 

(1987) 

, Tike 

car and Sin ah 


Komathiraj and Maharaj 


(2007). The 


spheroidal spacetime is found to accommodate superdense stars like neutron stars in both charged 


have been done by 


Sharma et al 

(2006), 

Sharma and Mukheriee 

(2001) 

Sharma anc 

Mukheriee 


(120021). Recently 


Pant and Maurya (2012) & 


charged f 


uid models have also been studied by 


Maurya et al (120151) . 


Maurya and Gupta (2011a,b,cj), 


In this paper, we have obtained a new class of solutions for charged fluid distribution on the 
background of pseudo spheroidal spacetime. Particular choices for radial pressure p r and electric 
field intensity E are taken so that the physical requirements and regularity conditions are not 
violated. The bounds for the geometric parameter K and the parameter a associated with charge, 
are determined using various physical requirements that are expected to satisfy in its region of 
validity. It is found that these models can accommodate a number of pulsars like 4U 1820-30, 
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PSR J1903+327, 4U 1608-52, Vela X-l, PSR J1614-2230, Cen X-3, given by iGangopadhyav e l al 


'( 20131) . When a _= 0, the model reduces to the uncharged anisotropic distribution given by 


Ratanpal et al ( 20151) . 


In section [2J we have solved the field equations and in section [3] we have obtained the 
bounds for different parameters using physical acceptability and regularity conditions. In section 
IH We have displayed a variety of pulsars in agreement wit the charged pseudo-spheroidal model 
developed. In particular we have studied a model for various physical conditions throughout the 
distribution and discussed the main results at the end of this section. 


2. Spacetime Metric 


We shall take the interior spacetime metric representing charged anisotropic matter 
distribution as 


ds 2 = e^dt 2 


l+K h 

i + 4 
^ w 


dr 2 — r 2 ( dO 2 + sin 2 9d(j ) 2 ), 


( 1 ) 


where K and R are geometric parameters and K > 1. T 


spacetime, has been s tudied by number of r e searc 


hers 


homas et al (2005); 


( 20ld : 


Thomas and Ratannall (120070 : 


lis spacetime, known as pseudo-spheroidal 


Tikekar and Thomas 


Paul et al (.20111 ): 


(1998 


1999 


2005b; 


Paul and Chattopadhyav 


Chattopadhyay and Paul ( 2012 ) have found that it can accommodate compact superdense 


stars. 


Since the metric potential g rr is chosen apriori, the other metric potential v (r) is to be 
determined by solving the Einstein-Maxwell field equations 

=s* (*;'+»>+£/)> <2) 

where, 

T/ = (p + p) iiiid - pdj, 


( 3 ) 
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7r/ = V3S CjC j - i [uiu j - <5/j , 


Ej = T + jF„F"5/) . (5 

Here p, p, Uj, S and c', respectively, denote the proper density, fluid pressure, unit-four velocity, 
magnitude of anisotropic tensor and a radial vector given by (o. —e A/2 .0. o). Fjj denotes the 


anti-symmetric electromagnetic field strength tensor defined by 


_ 

dxi d.Xj ’ 


which satisfies the Maxwell equations 


Fij,k + Fjk,i + FkiJ — 0 , 


F*y/=g)=47Cy/=gJ l , 


where g denotes the determinant of gij, A; = (0(r), 0,0,0) is four-potential and 


/' = au l , 


is the four-current vector where a denotes the charge density. 


The only non-vanishing components of F, j is Fo\ — —F\q. Here 

v+A 

Foi = — 6 9 [ Anr 2 oe X ^ 2 dr, 

F Jo 

and the total charge inside a radius r is given by 


q[r ) — An Or 2 e^/ 2 dr. 

Jo 


The electric field intensity E can be obtained from E 2 — —Fo\F m , which subsequently reduces to 


(12) 
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The field equations given by © are now equivalent to the following set of the non-linear ODE’s 


1 —e 


■ + ■ 


e~ k V 


= 8 np +E , 


(13) 


Ay' 


= 8 np r -E 2 , 


(14) 


-A 


where we have taken 


v" 


,'2 


2 + T" 


v'A' v'-A' 

- + ■ 


Pr = P + 


2 r 


2 S 

7 !’ 


= 87T/?j_ +£■“, 


(15) 


(16) 


Because e^ the metric potential A is known function of r. The set of equations (fl3l) 


i * K i 


P±=P 


S 


(17) 


1+ 


/? 

(fT5b are to be solved for five unknowns v, p, p r , and E. So we have two free variables for 
which suitable assumption can be made. We shall assume the following expressions for p r and E. 

K -1 1 


8 Kp r — 


r 

R 1 


R 2 


l + K W 


E 2 = CC(K~ 1 ) 


r 2 

PI 


R 2 


1+^ 


(18) 


(19) 


It can be noticed from equation (fl8l) that p r vanishes at r — R and hence we take the geometric 
parameter R as the radius of distribution. Further p r > 0 for all values of r in the range 0 <r < R. 
It can also be noted that E 2 is regular at r — 0. On substituting the values of p r and E 2 in (fl4l) we 
obtain, after a lengthy calculation 


e v - CR- 


\jP—( 2 +a)K+a+lj 


1 +K 


R 2 


( K±a+\\ 

f 2 K J 


K-a -3 


1 + 


R 2 


( 20 ) 
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where C is a constant of integration. Hence, the spacetime metric takes the explicit form 


ds 1 = CR 


\K 2 ~(2+a)K+a+\\ 


!K+a+l \ 
2 \{ 2K ) 


1+K W 


K-a -3 
\ - 2 - 


1 + ^2 I dt 


( 21 ) 


- -jr- J dr 2 — r 2 (dQ 2 + sin 2 dd<f> 2 ). 

1 + F / 

The constant of integration C can be evaluated by matching the interior spacetime metric with 
Riessner-Nordstrom metric 

-l 


ds 2 — l 1 — — + ) dt 2 — l 1 — — + \ ) dr 2 — r 2 {dO 2 + sin 2 0d(j) 2 ), (22) 


r r- 


2m q 


r r 


across the boundary r = R. This gives 


M = 


R [K 2 + a(K-l) -1] 


and 


C = R' 


2 (l+£)“ 

[+-(2+a)tf+a+l] /3Ar+a+l\ / a-K+5\ 

- K - (1 +A')~( 2K ) 2\ 3 ). 


(23) 


(24) 


Here M denotes the total mass of the charged anisotropic distribution. 


3. Physical Requirements and Bounds for Parameters 


The gradient of radial pressure is obtained from equation (fl8l) in the form 


8 = 
dr 


2r(K — 1) 1 + 2 K-K 


R 4 


< 0 . 


1+A- 


W 


(25) 


It can be noticed from equation (l25l) that the radial pressure is decreasing function of r. Now, 
equation (fl3l) gives the density of the distribution as 


&7tp 


K — 1 \ 3 + (^ — cc)^T 


( 26 ) 
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The conditon p(r — 0) > 0 is clearly satisfied and p(r — R) >0 gives the following inequality 
connecting a and K. 

0 < a <3+K. (27) 


Differentiating (l26l) with respect to r, we get 


dp 2r(K-l) 5K + a + K(K- a )^2 

8 71 —— — — 


dr 


R 4 


1+*F 


(28) 


It is observed that ^(r = 0) = 0 and c -^{r — R) <0 leads to the inequality 

K 2 - K(a - 5) + a >0. 

The inequality (l29l) together with the condition > 1 give a bound for a as 

K(K + 5) 


(29) 


0 < a < 


Z-l 


(30) 


The expression for is 


87Tp ± = 


4K -4+X l £+X 2 £+X 3 ^ 


R* 


(31) 


R 2 [4+y 1 ^+y 2 ^+y 3 ^+4k 2 ] 

where, X { = 4K 2 + (-l2oc-l6)K+\2a+l2,X 2 = 6K 3 + (-l0oc-22)K 2 + (4a+l4)K+6oc+2, 
X 3 = £ 4 + (-2a-4)£ 3 + (cr + 2a + 6)£ 2 + (-2a 2 -2a-4)£ + a 2 + 2a + l,yi = 12K + 4, 
y 2 = 12 K 2 + 12 K and Y 3 = 4K 3 + UK 2 . 

The condition p±> 0 at the boundary r — R imposes a restriction on K and a respectively 
given by 

K> 2^3 — 1 (32) 


and 


0 < a < 


10 + 5 K + K 2 IS9+102K + 51K 2 + SK 3 

('f-i) 2 ' 


(33) 


K- 1 
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The expression for is given by 

dpj_ ~ r (&K 2 + (12a + 8 )K — 12a — 16 +A\^ +A 2 ^4 4-A 4 ^^ 

dr R 4(2 + B l £_ + B 2 £ + B 3 ^ + B 4 ^ + B 5 j^ + 2K^') ’ ^ 

where, Aj = -4K 3 + (28 - 4a)K 2 + (16a - 20)K- 12a-4, A 2 = 3K 4 + (-4a - 4)K 3 + 
(-3a 2 - 28a - 30)£ 2 + (6a 2 + 44a + 36)K - 3a 2 - 12a - 5, A 3 = 10£ 4 + (-16a - 36)£ 3 + 
(-2a 2 + 4a + 16)£ 2 + (4a 2 + 16a + 12)^ - 2a 2 - 4a - 2, A 4 = K 5 + (-2a - 4)if 4 + 

(a 2 + 2a + 6)K 3 + (-2a 2 - 2a - 4)fc 2 + a 2 + 2a + 1, B x = 8k + 4, B 2 = UK 2 + \6K + 2, 

B 3 = 8 K 3 + 24K 2 + 8 K, 2K 4 + 16 K 3 + 12£ 2 and B 4 = 4K 4 + 8 K 3 . 


The value of - = 0 at the origin and d!2L (r — R) <0 gives the following bounds for K and 


a respectively 


2\/T3 — 5 < K <5 


(35) 


and 


K 3 + 10K 2 + 25K 

0 < a <-3- 

K 2 -6K + 5 


20 


■ + 


16 K 5 + 233 K 4 + 252K 3 + 218K 2 - 188K + 265 

I o 




(36) 


In order to examine the strong energy condition, we evaluate the expression p — p r — 2p ± at 
the centre and on the boundary of the star. It is found that 


(P — Pr~ 2p_i_) (r = 0) = 0, 

and (p — p r — 2p 1 _) (r — R) > 0 gives the bound on K and a, namely 

1 <K< 1+2V6 


(37) 


(38) 


0 < a < 


8 + 3 K + K 2 
K- 1 


+ 


/41+46K + 49K 2 + 8K 3 


( K-iy 


(39) 


The expressions for adiabatic sound speed 


dp r 

dp 




in the radial and transverse directions, 


respectively, are given by 

d p r 1 + 2 K — K^r 

dp 5k + a+K(K-a)£' 


( 40 ) 
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and 


( 1+jS: f) 

SK~ + (12a + 8 )K —12a—16 +Ci^t+C2^4+C3^+C4^ 


2 (^-l) 

1 1 
L/l 

+ 

R 

+ 

"R 

1 

R^ 

1_1 

2 + T>i^ 2 +Z>2^4 +Z>3^ +Z>4p + D$ JjjQ + 2K 4 

(41 


where, Ci = -4K 3 + (18 - 4a)K 2 + (16a - 20)K - 12a - 4, C 2 = 3K 4 + (-4a - 4)K 3 + 
(-3a 2 -28 a- 30) K 2 + ( 6a 2 + 44a + 36) K - 3a 2 - 12a-5, C 3 = 10 K 4 + (- 16a -36)A^ 3 + 
(—2a 2 + 4a +16) £ 2 + (4a 2 + 16a +12) £ - 2a 2 - 4a - 2, C 4 = K 5 + (-2a - 4)K 4 + 

(a 2 + 2a + 6) K 3 + (-2a 2 - 2a - 4) K 2 + (a 2 + 2a +1) K, D\ = 8^ + 4, D 2 = 12 K 2 +16tf + 2, 
D 3 = 8 K 3 + 24K 2 + 8 K, D 4 = 2 K 4 + \6K 3 + \2K 2 and D 5 = 4K 4 + 8 K 3 . 


The condition 0 < ^ < 1 is evidently satisfied at the centre whereas at the boundary it gives 


a restriction on a as 


0 < a < 


K 2 + 4K — 1 
K-\ 


, K> 1. 


(42) 


Further ^£4- < 1 at the centre will lead to the following inequalities 


K> 


(43) 


and 


0< a<-(3K-4). 

Moreover at the boundary (r = R ), we have the following restrictions on K and a. 


(44) 


—5 + 2\/T3 < K < 5 


(45) 


and 


0 < a < 


K 3 + 10K 2 + 25K-20 /16^r233^4^^3^^2r788Z+265 


K 2 — 6K +5 


- + 


{K 2 -6K + 5) : 


(46) 


The necessary condition for the model to represent a stable relativistic star is that T > ^ 
throughout the star. T > | at r — 0 gives a bound on a which is identical to (l27l) . Further, T —> °° 
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as r —>• R and hence the condition is automatically satisfied. It can be noticed that E = 0 at r — 0, 
showing the regularity of the charged distribution. 


The upper limits of a in the inequalities (l27l) . (l30l) . (1331) . (l36l) . (l39l) . (l42l) and (l44l) for different 
permissible values of K are shown in Tabled] It can be noticed that for 2.4641 < K < 3.7641 the 
bound for a is 0 < a < 0.6045. 
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Table 1: The upper limits of a for different permissible values of K. 


Inequality Numbers 


K 

m 


(133]) 

m 


m 

m 

2.4641 

5.4641 

12.5622 

0.0000 

10.1962 

1.6962 

0.0802 

30.9893 

2.5041 

5.5041 

12.4932 

0.0170 

10.1635 

1.7562 

0.0938 

30.6186 

2.6041 

5.6041 

12.3445 

0.0599 

10.0977 

1.9062 

0.1287 

29.7861 

2.7041 

5.7041 

12.2250 

0.1036 

10.0514 

2.0562 

0.1648 

29.0693 

2.8041 

5.8041 

12.1299 

0.1480 

10.0213 

2.2062 

0.2021 

28.4488 

2.9041 

5.9041 

12.0552 

0.1931 

10.0048 

2.3562 

0.2405 

27.9094 

3.0041 

6.0041 

11.9980 

0.2388 

10.0000 

2.5062 

0.2798 

27.4388 

3.1041 

6.1041 

11.9557 

0.2852 

10.0052 

2.6562 

0.3201 

21.0211 

3.2041 

6.2041 

11.9263 

0.3321 

10.0189 

2.8062 

0.3612 

26.6662 

3.3041 

6.3041 

11.9082 

0.3795 

10.0401 

2.9562 

0.4030 

26.3495 

3.4041 

6.4041 

11.8998 

0.4275 

10.0679 

3.1062 

0.4457 

26.0714 

3.5041 

6.5041 

11.9002 

0.4760 

10.1015 

3.2562 

0.4890 

25.8272 

3.6041 

6.6041 

11.9082 

0.5251 

10.1401 

3.4062 

0.5330 

25.6130 

3.7041 

6.7041 

11.9230 

0.5745 

10.1833 

3.5562 

0.5776 

25.4254 

3.7541 

6.7541 

11.9327 

0.5995 

10.2065 

3.6312 

0.6001 

25.3407 

3.7641 

6.7641 

11.9348 

0.6045 

10.2112 

3.6462 

0.6047 

25.3244 


4. Application to Compact Stars and Discussion 


In order to compare the charged anisotropic model on pseudo-spheroidal spacetime with 
observational data, we have considered the pulsar PSR J1614-2230 whose estimated mass and 
radius are 1.91 Mq and 9.69 km. On substituting these values in equation (l23l) we have obtained 
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the values of adjustable parameters K and a as K — 3.58524 and a = 0.292156 respectively which 
are well inside their permitted limits. Similarly assuming the estimated masses and radii of some 
well known pulsars like 4U 1820-30, PSR J1903+327, 4U 1608-52, Vela X-l, PSR J1614-2230, 
Cen X-3, we have displayed the values of the parameters K and a, the central density p c , surface 
density Pr, the compactification factor u — -jjp{r = 0) and charge Q inside the star in Tabled 


From the table it is clear that our model is in good 


data of pulsars given by 


Gangopadhyav et a l (2013). 


agreement with the most recent observational 


Table 2: Estimated physical values based on the observational data 


STAR 

K 

M 

R 

Pc 

Pr 


(tL 

Q 



(Ms) 

(Km) 

(MeV fm 3 ) 

(MeV fm 3 ) 



Coulomb 

4U 1820-30 

2.815 

1.58 

9.1 

1980.14 

250.46 

0.256 

0.461 

4.031 x 10 20 

PSR J1903+327 

2.880 

1.667 

9.438 

1906.90 

235.92 

0.261 

0.460 

4.184 x 10 20 

4U 1608-52 

3.122 

1.74 

9.31 

2212.22 

252.97 

0.276 

0.455 

4.127 x 10 20 

Vela X-l 

3.078 

1.77 

9.56 

2054.99 

238.25 

0.273 

0.456 

4.240 x 10 20 

PSR J1614-2230 

3.585 

1.97 

9.69 

2487.35 

248.17 

0.300 

0.448 

4.262 x 10 20 

Cen X-3 

2.589 

1.49 

9.178 

1705.08 

233.65 

0.239 

0.466 

4.044 x 10 20 


In order to examine the nature of physical quantities throughout the distribution, we 
have considered a particular star PSR J1614-2230, whose tabulated mass and radius are 
M — 1.97Mq, R — 9.69 km. Choosing K = 3.58524 and a = 0.292156, we have shown the 
variations of density and pressures in both the charged and uncharged cases in Figure [H Figure [2] 
and Figure [3l It can be noticed that the pressure is decreasing radially outward. The density in the 
uncharged case is always greater than the density in the charged case. Similarly the radial pressure 
p r and transverse pressure p± are decreasing radially outward. Similar to that of density, p r and 
p ± in the uncharged case accommodate more values compared to charged case. 
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The variation of anisotropy shown in Figure |4] is initially decreasing with negative values 
reaches a minimum and then increases. In this case also anisotropy takes lesser values in the 
charged case compared to uncharged case. The square of sound in the radial and transverse 
direction (i.e. ^ and ^-) are shown in Figure \5\ and Figure [ 6 ] respectively and found that they 
are less than 1. The graph of p — p,- — 2p± against radius is plotted Figure [7] It can be observed 
that it is non-negative for 0 < r < R and hence strong energy condition is satisfied throughout the 
star. 


A necessary condition for the exact solution to represent stable relativistic star is that the 
relativistic adiabatic index given by T — l>r ^ should be greater than The variation of 
adiabatic index throughout the star is shown in Figure [ 8 ] and it is found that T > 3 throughout 
the distribution both in charged and uncharged case. Though we have not assumed any equation 
of state in the explicit form p r — p r (p) and p± — we have shown the relation between 

PnP±_ against p in the graphical form as displayed in Figure [9] and Figure [TOl For a physically 
acceptable relativistic star the gravitational redshift must be positive and finite at the centre and 
on the boundary. Further it should be a decreasing function of r. Figure |TT] shows that this is 
indeed the case. Finally we have plotted the graph of E 2 against r which is displayed in Figure [I2l 
Initially E 2 increases from 0 and reaches a maximum values and then decreases radially outward. 


The model reduces to the uncharged anisotropic distribution given by 

a = 0 . 


Ratanpal et al ( 2015b when 
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Fig. 1.— Variation of density against radial variable r. 
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r(Km) 

Fig. 2.— Variation of radial pressures against radial variable r. 
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Fig. 3.— Variation of transverse pressures against radial variable r 
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Fig. 4.— Variation of anisotropies against radial variable r. 
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r{Km) 


Fig. 5.— Variation of against radial variable r. 
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Fig. 6.— Variation of against radial variable r. 
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Fig. 7.— Variation of strong energy condition against radial variable r. 



Fig. 8.— Variation of T against radial variable r. 
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Fig. 9.— Variation of pressures against density for charged case. 



Fig. 10.— Variation of pressures against density for uncharged case. 
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r (Km) 


Fig. 11.— Variation of gravitational redshift against radial variable r. 



r (Km) 


Fig. 12.— Variation of E 2 against radial variable r. 
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